Does the nuclear spin relaxation rate in superconductors depend on disorder? 
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We calculate the relaxation rate of a nuclear spin in s-wave superconductor with nonmagnetic 

■ impurities, including the strong-coupling effects. We show that in a weakly disordered three- 
' dimensional system the corrections due to disorder are negligibly small. 
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00 . I. INTRODUCTION 

(N 

In this article we revisit a long-standing problem about the effect of nonmagnetic impurities on the nuclear spin 
relaxation rate 1 R S (T) in isotropic superconductors. The first calculation of R s in a clean superconductor using 
Q ■ the Bardeen-Cooper-Schrieffcr (BCS) model was done by Hebel and Slichter 2 . They noticed that the expression 
O | for the relaxation rate is logarithmically divergent at all temperatures below T c , due to the singularity of the BCS 

■ quasiparticle density of states (DoS) at the gap edge. This singularity survives even in the presence of scalar disorder 
Qh' (since, according to Anderson's theorem, the superconducting gap is not affected by nonmagnetic impurities) and 

was removed in Ref^ by phenomenologically introducing some energy level broadening. A few years later, Maki and 
Fulde 3 calculated the relaxation rate in a superconductor with magnetic impurities, in which case the DoS singularity 
is smeared and the expression for R s is finite at all T. It was shown in Ref.— that the impurity vertex corrections to 
Ch ■ the relaxation rate are negligible. In a parallel development, the Hebel-Slichter arguments about the importance of 
I \ the level broadening were put on a quantitative footing by Fibich^, see also Ref£, who calculated the relaxation rate 
in the clean case using the Eliashberg theory of electron-phonon superconductivity^ (the so-called strong coupling 
model) and found that inelastic scattering indeed removes the gap-edge singularity and produces a finite relaxation 
rate. 

The effect of scalar impurities in the strong-coupling regime remained an open problem for another three decades 
until Choi and Melei presented a calculation of the relaxation rate, in which they questioned the accuracy of the Maki- 
Fulde result and claimed that the impurity vertex corrections are important. As a result, R s appreciably increases 
with nonmagnetic disorder, especially in the dirty case £ < £o, where I is the elastic mean free path and £o is the 
superconducting coherence length. In the present article we try to resolve this controversy. The article is organized 
\q ' as follows: In section [TTl we calculate the relaxation rate in the weak-coupling BCS limit, using the exact eigenstates 
method to perform the disorder averaging. In section Hill the exact eigenstates method is generalized to the strong- 
t-H 1 coupling case. In section HVl we calculate the relaxation rate in the strong-coupling case by another method, using a 
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^ II. EXACT EIGENSTATES METHOD: WEAK COUPLING 



■ Without loss of generality we consider the case of a nuclear spin 1=1/2 located at the origin of the crystal lattice. 
Higher values of / change only the overall prefactor in the expression for the relaxation ratei, which drops out of 
the ratio of the relaxation rates in the superconducting and the normal states. We assume, following Hebel and 
Slichter— that, while the resonance is observed in a strong field in the normal state, the relaxation takes place in 

, a uniform superconducting state after switching off the field. The spin-lattice relaxation rate due to the hyperfinc 

■ contact interaction of the nucleus with the band electrons is given by 

a; r^^=-j> hm ^ a) 

where J is the hyperfine coupling constant, uiq is the NMR frequency, and K is the retarded correlator of the 
transverse components of the electron spin density at the nuclear sitei^, which is calculated as follows. We introduce 
K(y m ), with v m = 2ixmT , as the Fourier transform of 

K(t) = -{{T T S+(0,T)S-(0,0))) imp . (2) 

Here (• ■ •) denotes the quantum statistical average and (■ • -)i mp denotes the averaging over impurity configurations 9 -^. 
The spin density operators in the Matsubara representation are given by S±(r, r) = e Hr S±(r)e~ Hr , where H is the 
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electron Hamiltonian, to be defined below, and 

S+W^}(r)fi(r), S_(r) = ^I(rM(r). (3) 

We use the units in which fi — kg — 1, and the spin quantization axis is chosen along the external magnetic field H . 
The retarded correlator in Eq. |TJ) is obtained by the analytic continuation: K r {uiq) = K(u m )\i Urn — >u> +iO- 

The properties of our system in the superconducting state can be described using a standard field-theoretical 
formalism in terms of the normal and anomalous Gor'kov functions 9 , which can be combined into a 2 x 2 matrix 
Green's function 

G(rr-r'r')-( G ^^ r '^') -F{r,T;r',r')\ 
1 ' ' ' ' ~ \ -F\r, r; r', r') -G(r>, r'; r,r) ) ' (4) 



Introducing the two-component Nambu operator: 



,ii 



*(r,r)= (JJ^j) , <t>(r,T) = (^(r lT ) ^(r,r)), (5) 

the matrix Green's function can be written in a compact form: 

G(r, r; r', r') = -(T T *(r, r)*(r', r')). (6) 

The Matsubara spin correlator ^ can be expressed as the impurity average of the product of two matrix Green's 
functions: 

Hvm) = \ T Y1 Tr(G , (0,0;c; n + ^ m )G(0,0; W „)) lmp , (7) 

n 

where lu u = (2n + 1)ttT. Note that this representation of the spin correlator in terms of the G's is not unique, see 
section IIVI below. 

We consider a single spin-degenerate electron band in a three-dimensional (3D) crystal and neglect the spin-orbit 
coupling. The electron Hamiltonian is written as H — Hq + Hi nt , where 

H = [ d 3 r ^ a (r)hoMr) (8) 



describes noninteracting electrons, with a =|, I being the spin projection, and 

ho - -P- + U(r) - e F . (9) 

Without loss of generality we assume an isotropic parabolic band. The impurity potential U(r) is characterized by 
the correlator 

(U( ri )U(r 2 )) mp = ^-l— 5( ri - r 2 ), (10) 

where Np is the density of states at the Fermi level, and r e ; is the elastic mean free time. 

In a weak-coupling BCS superconductor, the interaction Hamiltonian which describes s-wave singlet pairing has 
the following form: 

H int = \J d 3 rV a ^ jS iPi(r)^l(r)^(r)Mr), (11) 

where V a p n s = — (A/2)(i<T2) a; 3(i(T2)^ 5 , and A > is the coupling constant. Treating the pairing interaction in the 
mean-field approximation, the equation of motion for the matrix Green's function (|4|) can be written in the form 

~%r) -t+l ) 6(r ' T; r ' T ' } = foS (r - r ' )S{T - /} ' (12) 
where fo is the unity matrix in the Nambu space, and A(r) = XF(r, r; r, r) is the gap function. 
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The spin correlator will now be calculated using the exact eigenstates method, see e.g. RefJ^, which allows one 
to relate properties of the system in the superconducting state to those in the normal state. The key assumption is 
that the disorder-induced fluctuations of the gap function can be ignored: 

A(r) -> A = (A(r)) imp . (13) 

According to Anderson's theorem, the average gap is not renormalized by nonmagnetic impurities, i.e. the value of 
A at any temperature is the same as in the clean case. We introduce the exact eigenstates and eigenvalues of the 
single-particle Hamiltonian with impurities: 

hoVa(r) = Ca^a(r), J d 3 r \<p a (r)\ 2 = 1, (14) 

where, according to the definition ([9]), the chemical potential is included in £ a . Then the Gor'kov equations (fT2|) can 
be solved, with the result 

G{r,r';uj n ) = ^ g a {uj n )ip a {r)(p* a {r'), (15) 

a 

where 

5qK) = «a+g + A» ' (16) 

Substituting this in Eq. (|7|) we obtain, after the Matsubara frequency summation and the analytic continuation 

Im^M = lM0)l>bm 2 6(Ea _ Eb) 

OJ ^0 LJQ 2 \ L ' 

\ a.b 

SfU l+ f^)\ , (17) 



dE a I \ E a Eb 



where f(x) = (e x ^ T + 1) 1 is the Fermi function, and E a = yfQ + A 2 . 
Next, we introduce the local DoS in the normal state: 

tf(r,e) = X> (r)| 2 *(e-£ ), (18) 

a 

perform the disorder averaging in Eq. (|17[) . and insert the result in the expression ([T]) for the relaxation rate: 

R=J 2 N 2 F J de 1 de 2 T{e 1 ,e 2 ;u;)TZ{e 1 -e 2 ). (19) 

The superconducting gap enters here only through the function T, which in the weak-coupling case has the form 
T = T wc , with 

7T / £l£2 + A 2 ^ V^f+A 2 

r w ( ei> « ;w ) = 2 i + ^r ^ 



: [*(ei - e 2 ) + 8(e 1 + e 2 )]5(w -yje\ + A 2 ), (20) 



while all effects of impurities are contained in the local DoS correlator 

1 



n(n) = - T (N(0,e + n)N(0,e)) mp . (21) 



N 2 



In the normal state, we set A = in Eq. (|20|) and obtain R n = ir J 2 NplZ(0). To calculate the local DoS correlator, 
we write N(0, e) = -[G R {0, 0; e) - G A (0, 0; e)]/2m, where 

g^)(0,0;e)=^ (22) 
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are the exact retarded (advanced) Green's functions of the normal metal. Neglecting the difference between the local 
and the Fermi-level densities of states, we replace (N(0, e)),- mp = Np, and obtain 

K(Q) = 1 +STZ(n), (23) 

where 



5n(Sl) = ^2 j^t Re i(G R (0, 0; e + fl)G A (0, 0; e)) lmp 



-{G H (0, 0; e + ft)) lmp (G A (0, 0; e)) lmp ] (24) 

represents the impurity vertex corrections to the product of the retarded and advanced Green's functions. The terms 
containing the averages of two retarded or two advanced Green's functions vanish. The subsequent steps are standard. 
Introducing the disorder-averaged Green's functions 

G R W(k,e)= ] , (25) 

where = k 2 /2m — ep, the impurity ladder diagrams corresponding to S1Z can be summed, with the following result: 

^(0) = ^-Re f£± * { «;V (26) 
v ; irN F J (2tt) 3 1 - $(q, fi) ' v ; 

where 

1 f d 3 k R ( q 



■*^ic;/H a, ( l+ H«'(»-j 



l 



1 - iflTel + iT e iv F (k)q 



(27) 



vp(k) is the Fermi velocity, and the angular brackets denote the Fermi-surface averaging. At small il and q, &(q, fi) ~ 
1 + iflT e i — Dq 2 r e i, where D = v F T e i/3 is the diffusion coefficient. Therefore the impurity vertex corrections lead to 
a diffusion pole in the momentum integral on the right-hand side of Eq. (|2l)|) . In 3D, due to the convergence of the 
integral at q — > 0, the dependence of S1Z on f2 is not singular and can be neglected. On the other hand, <&(<?, fi) decays 
slowly at q — > oo: 

i 1 - ifir e i + iv F qT e i g-»oo 7r 1 

$(<7, S2) = In ~ , 

2vFqr e i 1 - ii2r e ; - ivFqr e i 2vft e i q 

for a 3D spherical Fermi surface. Thus it is necessary to introduce the ultraviolet cutoff of the order of the Fermi 
momentum hp in the integral (|26p . which gives the following estimate: 



-W?<<>) x £L q 2 dq J—l- oc (2S) 

Jo w F T e1 ? e ^ T e' 



and 

K(S1) = 1 + O ( -—) . (29) 

Although there is some enhancement of the local DoS correlator and, therefore, of the nuclear spin relaxation rate 
due to the diffusive motion of electrons, the magnitude of this effect in a weakly disordered (eFi~ e i 3> 1) 3D metal 
turns out to be negligibly small. Similar conclusions have also been reached in Ref. 13 . Thus we finally arrive at the 
following result in the normal state: 

R n =vJ 2 N 2 F , (30) 

which is known as the Korringa lawi. 

In the superconducting state, the integral of T wc on the right-hand side of Eq. (fl"9]) has a logarithmic divergency, 
whose origin can be traced to the square-root singularity at the gap edge in the BCS density of states. This singularity 
is smeared and the divergency is removed in the presence of magnetic impurities or if the gap is anisotropic 14 . In this 
work we focus on the gap smearing due to the strong-coupling effects in the Eliashberg formalism. 
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III. EXACT EIGENSTATES METHOD: STRONG COUPLING 



In this section, the exact eigenstates method is generalized to include the effects of electron-phonon interaction, see 
also Refs i 15 i 16 . The Hamiltonian is written as H = H Q + H p h + H e _ p h, where H describes noninteracting electrons 
and is given by Eq. {SJ, 

H ph = J2^(q)bljb qj (31) 

is the Hamiltonian of free phonons, with j labelling the phonon branches and u>j{q) being the phonon dispersion 
(recall that h = 1 in our units), and 

H e - P h= 9j(ki,k 2 )(b qj +b t _ qj )c{ il7 c k2(y , (32) 

ki,k 2 ,j 

describes the electron-phonon interaction (q = ki — k 2 ). The electron-phonon vertex has the following properties: 
<7j(fei, k 2 ) = g*(k 2 -, k\) due to the hermiticity of H e - p h, and gj{k\, k 2 ) — g*(—ki, —k 2 ) due to time-reversal symmetry. 
We assume unit volume and neglect the effects of impurities on the phonon spectrum. We also neglect the disorder 
effect on the electron-phonon vertices, which can be justified for long-wavelength acoustic phonons by the fact that 
electrons move with the lattice, including the impurity atoms, in order to preserve charge neutralit y 17 ' 18 . Although 
this argument cannot be extended to short- wavelength phonons, we make the usual assumption that the gj's are 
nonrandom functions. 

The matrix Green's function which includes the electron-phonon interaction but not disorder averaging satisfies the 
following equation: 

G = Go + Go^phG. (33) 

Here Go is the matrix Green's function of the disordered normal metal without phonons, whose coordinate represen- 
tation can be obtained by setting A = in Eq. (fT5|) : 

iuJ nTQ +£qf 3 



G (r,r ;u n ) = - ^ 2 <p a (r)<p a (r ), (34) 



and T, p h is the self-energy due to the electron-phonon interaction: 

E P h(k,k';cj n ) = -T}^ 5 kl ^w v k-k'Dj(k-ki,uj n -u; n ') 

n' ki,k' t ,j 

xgj(k, fei)5 J (A;i,fc')' f 3G(fci,fci;w„0f3, (35) 

where 



",■<</■"->' Tpr§^ ^'^ 



is the bare phonon propagator. The contributions to the self-energy from the diagrams with crossed phonon lines are 
neglected based on Migdal's theorem. 

The self-energy (|3"5"j) is still a random quantity. To make progress one has to assume that it can be averaged 
independently from the Green's functions, which amounts to replacing the random self-energy in the equation (|33|) 
by its disorder average: 

T, ph (k,k';oj n ) -> Ti ph (k,u) n ) = (T, ph (k,k';u} n )) imp . (37) 

This approximation is justified by the slow spatial variation of the self-energy compared to that of the electron Green's 
functions^. We use the subscript "ph" to emphasize that Y* p h represents only the phononic part of the full self-energy 
(the latter contains also the impurity part: S = T, p h + Si mp ). The Coulomb interaction can be included in a similar 
fashion. Although there are impurity vertex corrections to both the electron-phonon and the Coulomb vertices due 
to the diffusive motion of electrons, in a weakly disordered 3D superconductor those are smal l 15 ' 19 . Thus, the average 
self-energy can be written as 



1 r°° 

m F k' ^° 



■m 



f n 2 

XT 3 G(k',u) n - v m )f 3 , (38) 
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where 

a 2 F(fc,fc';ft) = N F ^\ gj {k,k')\ 2 5[VL - uj 3 {k - k 1 )], (39) 

3 

and G(k,u) n ) is the disorder-averaged Green's function of electrons. 

For the exact eigenstates method to work we have to neglect the anisotropy of the electron-phonon interaction, 
which gives T, P h(k, u) n ) = £ p /j(cJ n ). Similar to the usual decomposition of the full self-energy, £ — iu! n (l — Z)tq + 4>ti— , 
we represent its phononic part in the form 

S p fc(w„) = iaj n [l - Z ph (oJ n )]to + (j) p h{u n )Ti, (40) 

where Z p % and 4>ph are real and even functions of the Matsubara frequency. By analogy with the gap function 
A(oj n ) ee (f)(uj n )/Z(uj n ), one can define A ph (u n ) = c/) ph (uj n )/Z ph (uj n ). 

Replacing E p ^(fe, k! ', uj n ) — > S p (,(w„) in Eq. (|33p we find that the electron Green's function before disorder averaging 
has the form (|15|) . where g a is now given by the following expression: 

9a{u n ) = 3-=2~? N i C2 i 12 ? \ ' ( 41 ) 

Inserting this in Eq. Q , summing over the Matsubara frequencies and averaging with respect to disorder one obtains 
the expression (fT9|) for the relaxation rate, with T = T sc , where 

7T 

r sc (ei,e 2 ;o;) = - Tr [p(e x , w)p(e 2 , w)J, (42) 



„, v 1 T iw n Z„ft(o; n )ft) + ef 3 + <^ p /i(w n )fi 
p{e,uj) = im TTPn — % o pt - ? — ; — 

Using Eq. (|29|) . the integrals over ei and £2 can be calculated separately: 



(43) 



1 f°° ^ iuj n Z ph (uj n )T + ef 3 + p ,,(a;„)Ti _ iw„f + A(w„)fi 



TT ./_«, w'Z^Wn) + e 2 + </> 2 ,>„) ^2 + A2( Wn ) 

Here we used the fact that the gap function is not renormalized by impurities: A p h(uj n ) = A(uj n )2Q. Now one can 
perform the analytic continuation and obtain: 

d e p{ e,.) = -Re»;° + A ^\ (44) 

where the branch of the square root is chosen such that its real part has the same sign as u>, and A (a;) is the complex 
gap function for u> just above the real frequency axis 2 ^. Finally, using the normal-state relaxation rate (|30p. we obtain: 

Rn J-00 \ ^ 



Re - ) +(Re- AM 



y/u 2 - A 2 ( W )/ V y / LJ 2_ A 2( w ) 



(45) 



which coincides with the clean-limit expression derived by Fibich in RefA This shows the absence of the impurity 
effects on the nuclear spin relaxation rate in a weakly disordered strong-coupling superconductor [up to the terms of 
the order of (ei?r e z) -1 ], which is the main result of this article. 

Note that in general the disorder enters the expression for the relaxation rate through both the local DoS correlator 
1Z and r sc , see Eq. (|42|) . The latter depends on Z p h and <p p h, which are both renormalized by disorder. It is the 
smallness of nontrivial disorder-induced correlations of the local DoS in 3D, i.e. the fact that lZ(ei — 62) — 1, that 
allows one to express the relaxation rate entirely in terms of the gap function A(w), which is not affected by disorder. 
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IV. PERTURBATION THEORY IN THE LADDER APPROXIMATION 



As an additional check of our result (|45|) . we now calculate the spin correlator ([2]) using a direct summation of the 
impurity vertex corrections in a standard diagram technique. While some of our intermediate results look similar to 
those of Refii, the final conclusion turns out to be qualitatively very different. We find, in agreement with the exact 
eigenstates method of section IIII1 that the leading order correction to the nuclear spin relaxation rate resulting from 
impurity scattering is of the order of (cpT^)" 1 in bulk superconductors. 

We use the two-component Nambu operators ([5]) in the momentum representation: 

Cfe ^ = (cL fe T( ("r)) ' ^( T ) = (4 T W C -H(r))- (46) 
Then the spin-density operators ([3]) take the following form 

S+(0,t) = - Y / Ckr{r)(r + ) rs C k ' s (T), 
fe.fc' 

5 -(°>°) = 5$3Ck-(0)(r-)r.Cfc/ a (0) J ( 4? ) 
fc,fc' 

where r, s — 1,2 are the Nambu indices, and f± = f± ± if 2 . The Hamiltonian of the system contains the electron- 
phonon and the screened Coulomb interactions as well as the term describing the nonmagnetic impurity scattering. By 
applying the Wick theorem to the correlator ([2]), in which S± are given by the expressions (|47p . one finds in the clean 
case K(t) = (1/2) Ylkk' ^ r — r)(if2)G(fc', — r)(ir2)]. In the presence of impurities we obtain the following form 
for the disorder-averaged K(v m \ which includes only the second Born approximation impurity ladder diagram a 9 i 10 : 

-K + v m )) 



2 

k.q 



xt(q,u n ,u n + v m )], (48) 



where G is the matrix Green's function of electrons, disorder-averaged and fully dressed by all interactions. We assume 
a standard isotropic strong-coupling superconductor, for which 



i u n Z(w n )f + gfef 3 + (j>{u!n)T l 
U J lZ*{LO n )+e k + <f{Un) 



u(K,u) n ) — ..ir/n,. \ , e2 i ±11. . \ > v 4y / 



with the momentum-independent renormalization function Z(u> n ) and the pairing self-energy <p(u) n ), which are even 
functions of uj n . Both Z and contain the effects of impurity scattering, which drop out of the gap function 
A(u> n ) = (f>(u> n ) / Z (u> n ) in the isotropic single-band case^. The vertex function T satisfies the equation 

f (q, cu n , uj n + v m ) = it 2 + — — \ V" f 3 G(k + q, —w n - v m ) 

k 

xT(q,uj n ,uj n + v m )G(k,uj n )f 3 . (50) 

Note that Eq. (|48|) could also be obtained from the spin correlator ([7|), using the identity (if2)G(0, 0, — <J n )(iT>j) = 
G(0,0,cj n ). In the clean limit we recover from Eqs. (|48|) and (f50|) the Fibich's expression^ for the nuclear spin 
relaxation rate in isotropic single-band strong-coupling superconductors. 

In order to solve Eq. (|50|) . one represents the vertex function as a linear combination of the Pauli matrices in the 



Gor'kov-Nambu space: T = ^~\._ Tifi, where I\ = Ti(q, ui n , oj n + v m ) = (1/2) Tr (fir), and finds a set of four coupled 
algebraic equations: 

3 

r i = - ^2 Li J r J 

V- ( 51 ) 

r 2 = % - 2^ l 2jTj 

3 



r 3 = ^i 3j r j 
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where L{j = Lij(q, u n , u> n + v m ), with i, j = 0, 1, 2, 3, are defined by 



L 1 1 

lJ 2^r e; 2 



- Tr [fiG(fe + q, -w n - u m )fjG(k, u n )\. 



(52) 



One can see that L w = Lot, L 20 = -L a2 , L 21 = -L 12 , L 3Q = L Q3 , L 31 = L 13 , L 32 = -L 23 , so that only ten out of 
sixteen L^'s in Eq. (IFl]) have to be computed. The sum over k in Eq. (f5"2")l is calculated approximately by using 

iy !■! r+oo 

/ ds / <*£i>(6f + Q(«))> 

assuming that p is decreasing fast enough with Here Q(s) = Ae F {q/2k F ){q/2k F + s). In this way we find 



-^oo — f- 



/00' + $$' 
V DD 7 

/oo' + $$' 



1 , ioi = */- 



- 0'$ 



DD 



ill = /- 



V DD' 



- 1 



L 13 = / . [ - - — 



-^23 — /- 



D D 

o$' + o'$ 

DD' 



<>' 



(53) 



V DD' 
,'00' - $$' 

^33 = /- + 1 



where 



and 



n = uj n Z{uj n ), $ = <f>(u n ), L'=v / J + ¥2, 

O' = (ui n + v rn )Z{w n + u m ), $' = (f)(ui n + v m ), D' = ^Q' 2 + <P', 



/+ 
f- 



1 l^x 2 (x + l) 2 + {D + D') 2 /16e 2 F 



•In- 



32e F T e ix x 2 (x-l) 2 + (D + D') 2 /16ej 



1 1 



16£FT e ; X 



tan 



D + D' 



— tan 



D + D' 



&e F x(x — l) Ae F x(x + l) 

with a; = q/2k F . We note that the relationship between our /± and (f±) of Refii is (/+) = 8if + , (/_) = — 8/_, with 
the additional difference that our /± contain (D + D')/Ae F , instead of (D + D')/e F . 

We solved the system (I5T1) with the Ly's defined by the expressions |53|) , using MAPLE. The result is 



r„ 

Ti = - 
r 2 = i 



(l-2/_)2+4/2 V-D ' D> 



0' 



(1 - 2f-) 2 +4f 2 \D D' 

f--2(fl + fl) / OO'-tjKfr' 

(l-2/_) 2 + 4/ 2 ^ DD' 
. /_ -2(/ 2 +/2) ft$' + 0'$ 



(54) 



(I-2/_) 2 + 4/2 DD' ' 
Inserting these into the vertex function in Eq. (|48|) . we finally obtain: 

oo 

K{ Vrn ) = -T J2 A nBn, 
n— — oo 

w„(w„ + z/ m ) - A(uj n )A(uj n + v m ) 



(55) 



where 



yju 2 + A 2 (w„)a/(w„ + V m ) 2 + A 2 (iJ n + V m ) 



B n = 2nT e iN F 



f--2 



n 



/-(/--2(/| + / 2 )) 

(l-2/_)2+4/| + (l-2/_) 2 + 4/ 2 
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Clearly, in the limit r, 



oo only the first term in the square bracket in B n survives, and it is possible to integrate 

p 2 )/2, where p = (D + D')/ep. For all the terms in the sum 
55)) for which A n is nonzero [note that A(w„) ~ for \ui n \ greater than 10 times the maximum phonon 



it over q analytically, with the result Tr 2 N F ^J~(l 
over n in Eq 



bat L 

frequency], p can be set equal to zero. In this way one recovers Fibich's formula^ for the relaxation rate, after the 
sum over n is performed, followed by the analytic continuation iv m —> ujq + zO, in the limit luq — * 0. We note that our 
B ni containing the impurity vertex corrections, is different from the expression for the vertex corrections obtained by 
Choi and Mele 7 -, which we have not been able to reproduce. An additional difference is that we calculate analytically 
the momentum integrals which were treated in Ref. 7 using some approximation. 

Even for a finite r e ; one can replace p — > 0, which makes it possible to integrate over q the second and third terms 
in the expression for B n . Our final result for the nuclear spin relaxation rate in the superconducting state to the 
leading order in (eF^e;) 1 1S 



R, = 



nJ 2 Nl 



Re- 



1 



3tt 1 

16 £ F T e l 
UJ \2 



OUJ 



y^ 2 - A 2 (w) 



Re 



A( w ) 



V^ 2 -A 2 H 



(56) 



We see that the impurity vertex corrections turn out to be of the same order as the diagrams with crossed impurity 
lines, which we have neglected, and the ratio R s /R n is therefore unaffected by impurity scattering at rj 1 <§; €f- This 



condition is much weaker than I 3> £o, he. r e 
Refi 



<C Aq (Aq is the superconducting gap at zero temperature), found in 



V. CONCLUSIONS 



We conclude that the answer to the question in the title is negative. Using two different techniques, the exact 
eigenstates method and the usual diagrammatic perturbation theory in the ladder approximation, we have shown 
that the contribution of nonmagnetic impurities to the nuclear spin relaxation rate in a bulk superconductor with 
isotropic pairing is of the order of (ei?r e ;) _1 , i.e. very small. 
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